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We study the phenomenology of supersymmetric models in which gauge-singlet 
scalars mix with the MSSM sneutrinos through weak-scale A terms. After reviewing 
the constraints on mixed-sneutrino dark matter from measurements of Q*cdm an d 
from direct-detection experiments, we explore mixed-sneutrino signatures relevant to 
the LHC. For a mixed-sneutrino LSP and a right-handed slepton NLSP, decays of the 
lightest neturalino can produce opposite-sign, same-flavor (OSSF) dileptons with an 
invariant-mass distribution shifted away from the kinematic endpoint. In different 
parameter regions, the charginos and neutralinos produced in cascades all decay 
dominantly to the lighter sneutrinos, leading to a kinematic edge in the jet-lepton 
invariant-mass distribution from the decay chain q — > x~Q ~> &*lq, without an OSSF 
dilepton signature. We explore the possibility of using mass estimation methods to 
distinguish this mixed-sneutrino jet-lepton signature from an MSSM one. Finally, we 
consider signatures associated with Higgs-lepton or Z-lepton production in cascades 
involving the heavier sneutrinos. 

I. INTRODUCTION 

The overwhelming evidence that the majority of the matter in the universe is non-baryonic 
compels us to consider extensions of the standard model which include new fields that are 
electrically neutral. As was pointed out by Goodman and Witten [TJ fields at the weak 
scale naturally yield the appropriate relic abundance to explain the observed density of dark 
matter in the universe. Taken alone, this may be the strongest motivation for new physics 
at the weak scale. 
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A second element of physics beyond the standard model that has become firmly es- 
tablished in recent years is neutrino mass. Although it is usually assumed that neutrino 
mass is generated at very short distances, thus explaining its smallness through the seesaw 
mechanism (2j [3j, this hypothesis remains untested, motivating us to consider alternative 
possibilities. 

For example, in supersymmetric theories with right-handed neutrinos at or below the 
weak scale, small Dirac neutrino masses can be generated as a supersymmetry-breaking 
effect, or else small Majorana neutrino masses can be generated radiatively or through a 
weak-scale seesaw [U El E]- Moreover, if the scalar partners of the right-handed neutrinos mix 
appreciably with the MSSM sneutrinos, the lightest "mixed" sneutrino can be considerably 
lighter than the Z boson [U EJ [7j, and can also be a viable dark matter candidate jl]. 

The existence of new states beyond those of the MSSM may be crucial for searches at the 
LHC. Unlike in previous experiments, new particles are likely to be produced in bunches in 
potentially long cascades. The presence of new states in the cascade chains can lead to new 
signatures and remove expected ones. It is an intriguing and important question to what 
extent analyses that can be performed at the LHC might give some indication that a mixed 
sneutrino is present in the spectrum of the theory. 

In this paper we study the cosmology and potential LHC signatures associated with 
mixed sneutrinos. In section [II], we consider the possibility of mixed-sneutrino dark mat- 
ter, calculating the relic abundance for a range in parameters, and imposing all present 
constraints from direct detection experiments, including the most recent from XENON [8]. 
In agreement with [9], we find that mixed-sneutrino dark matter is viable over substantial 
parameters regions. We consider the extent to which lepton number violation in the sneu- 
trino mass matrix might suppress rates at direct detection experiments, and discuss the 
connection to neutrino masses. 



In section III we explore the collider phenomenology of mixed-sneutrinos. The first sig- 
nature we consider, which persists even for very small sneutrino mixing angles, arises from 
leptonic decays of the lightest neutralino in the case where the LSP is the lightest sneu- 
trino and the NLSP is a right-handed slepton. The possibility of lepton production from 
the decays of the lightest neutralino was also pointed out in ref. [10], which studied the 
phenomenology of a sneutrino NLSP with a gravitino LSP. We study the invariant- mass 
distribution of opposite-sign, same flavor dileptons from these decays and find that it is 



3 



shifted away from the kinematic endpoint. The other signatures we consider require some- 
what larger mixing angles, as they involve decays of non-NLSP superpartners straight to 
the lightest sneutrino. For example, in a broad region of parameter space, the gauginos pro- 
duced in cascades decay almost exclusively directly to the lightest sneutrinos. In this case 
one has a kinematic edge in the lepton-jet invariant mass distribution, without a correspond- 
ing dilepton edge. We discuss the possibility of distinguishing this signature from MSSM 
ones, for example by using recently proposed methods to estimate the masses involved in 
cascade decays. Finally, we consider signatures from Higgs-lepton or Z-lepton production in 
cascades involving the heavier sneutrinos, which can lead to distinctive bbl, 77Z or trilepton 
invariant mass distributions. 

II. MIXED SNEUTRINO DARK MATTER 

To be viable, WIMP dark matter candidates must pass three essential tests. First, they 
must be neutral, both to allow early growth of structure and to have evaded detection. 
Second, their relic abundance must match the measured value of the dark matter energy 
density, VL CDM h 2 ~ 0.1 [TTj . Third, given the appropriate relic density they must evade 
direct-detection experimental limits, the most severe of which presently come from XENON 
[E] and CDMS [TJ. 

The sneutrino was long ago considered an intriguing dark matter candidate [131 E] > but is 
no longer viable as it fails the combined relic abundance and direct detection requirements. 
In particular, a light sneutrino with the appropriate relic abundance would significantly 
modify the invisible Z-width, in conflict with observation. A heavy sneutrino must be of the 
order 600 GeV to achieve the correct relic abundance [T5], and even then is in clear conflict 
with direct detection experiments. Similarly, even if a moderate-mass (~ 100 GeV) sneutrino 
had come out with the correct relic abundance, it would have been seen at experiments 
such as CDMS and XENON. One proposal for saving sneutrino dark matter is to suppress 
coannihilation of the sneutrino's scalar and pseudo-scalar components by making them non- 
degenerate [in], which also eliminates direct detection constraints arising from Z-exchange 
contributions to the scattering of sneutrinos off of nuclei. Unfortunately, this scenario implies 
a v T mass well above the experimental limit. 

The problems with sneutrino dark matter mainly stem from the large coupling of sneutri- 
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nos to Z bosons. However, because the sneutrino is neutral under electromagnetism, it is free 
to mix with any additional neutral scalar field, assuming that the field carries lepton number 
or that lepton number is not a good symmetry of the low-energy theory. This possibility was 
explored in jU El [7]. The mixing suppresses the coupling of the lightest sneutrino to the Z, 
and its mass is allowed to be less than for mixing angles satisfying sin8 J$ 0.4. Because 

the sneutrino annihilation rate in the early universe is also suppressed, the appropriate relic 
abundance can be achieved [I]. 

Related scenarios for sneutrino dark matter include non-thermal right-handed-sneutrino 
dark matter (where the mixing is extremely tiny) [TTJ [TS], and thermally produced right- 
handed-sneutrino dark matter in the presence of an extra U(l) [19|. 

The outline for the rest of the section is as follows. First, we review models of mixed- 
sneutrino dark matter. Then we discuss the relic abundance calculation and identify cos- 
mologically preferred parameter regions. With these results in mind, we review constraints 
from direct-detection experiments, and find, in agreement with [5], that significant regions 
of parameter space remain viable. Lepton-number violation in the sneutrino mass matrix 
can suppress the scattering of sneutrinos off of nuclei, and thus direct-detection rates, but it 
also radiatively generates neutrino masses that tend to be beyond experimental limits. We 
discuss a few scenarios in which these neutrino masses are not problematic, and then briefly 
consider the implications of sneutrino dark-matter for neutrino telescope indirect detection 
experiments. Finally, we comment on scenarios in which the gravitino is the LSP, with a 
mixed-sneutrino NLSP. 



The model we consider is quite simple. To the MSSM, we add one or more additional 
standard- model-singlet superfields iVj, with supersymmetry-breaking trilinear couplings of 
the form A^hil^h u . Restricting ourselves for the moment to one generation, this leads to a 
mass matrix of the form 



with mass eigenstates V\ = cos On* — sin#z> and v<i = sin#7i* + cos 6v. Motivated by the 
possibility of mixed-sneutrino dark matter, we take sin 2 # < 0.5, so the lighter state is 



A. Mixed Sneutrinos with Large or Small Yukawas 




(1) 
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more singlet than active sneutrino. If this lighter state is heavier than rriz/2 there are 
no immediate constraints on sin#, while if it is lighter than mz/2, the Z-width constraint 
requires sin 9 < 0.4. 

The A terms for the superpartners of the standard model fermions are typically thought 
to be related to the associated Yukawa couplings. Given the apparent smallness of the 
neutrino Yukawa couplings, one thus might not expect sizeable mixing between the active 
and sterile sneutrinos. However, it is possible that bare Yukawa couplings for the neutrinos 
are forbidden by aZ7(l) n ®i7(l)j symmetry that acts independently on the singlet and lepton- 
doublet superfields. If this symmetry is broken only by supersymmetry-breaking fields, then 
weak-scale A terms and tiny Yukawa couplings are perfectly compatible. 

Moreover, one can instead have large Yukawa couplings and still have massless neutrinos, 
as we now describe. If the fields Ni come with fields Ni which carry opposite lepton number 
charge, we can consider the following superpotential, 

W D XNLH U + m N NN. (2) 

When the Higgs acquires an expectation value, there are Dirac masses between v and n, 
as well as between n and n. Because of the mismatch between states with lepton number 
+1 and —1, there is a massless state in the theory. This is essentially the same mechanism 
that keeps the neutrino light in the standard model. If > Xv u the massless state will 
then be mostly standard model-neutrino. Constraints on this scenario come from a variety 
of precision electroweak measurements, principally from measurements of the couplings of 
charged leptons to neutrinos. For light (m^ ~ rnz) neutrinos, Xv should be smaller than 
about m T . However, for heavier neutrinos, a larger Yukawa is allowed, even for much lighter 
sneutrinos. This setup results in a 3 x 3 sneutrino mass matrix instead of the 2x2 one of 
eqn. (TO), 

m 2 L + \m 2 z cos 2/3 Av sin (3 + A/if cot (3 Xtunv tan j3 
Av sin (3 + A/if cot (3 m| • (3) 

^ Xm N vtan(3 m| j 

In the limit where (and thus m|) are very large, the effective mass matrix for the 
lighter sneutrinos is the same as in eqn. (JI]) except with the replacement A — > X, where 
X = A + X/2cot(3. 

In this setup, non-zero neutrino mass can be generated through higher-dimension opera- 
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tors, or radiatively if a small lepton-number- violating terms appear in the full sneutrino mass 
matrix. For our calculations of direct detection rates and relic abundances we will restrict 
ourselves to the model with negligible Yukawa couplings, and leave a thorough analysis of 
the relic abundance of the Yukawa model to future work. However, it is worth noting that 
even within the mixed-sneutrino framework, there is great room for variation. 

B. Relic Abundance of Mixed Sneutrinos 

The dominant annihilation channels for mixed sneutrinos in the early universe are shown 
in figure [I] These include s-channel Z exchange, t-channel neutralino exchange (to vv or 
vv), and s-channel Higgs exchange (to fermions, or, for heavier sneutrinos, to gauge bosons 
and Higgs bosons). The contribution from Higgs exchange is enhanced by the large A-terms, 
and is often dominant. 



V V 




FIG. 1: Diagrams contributing to the mixed-sneutrino annihilation rate in the early universe. 

To calculate the relic abundance of the mixed-sneutrino LSP, we use the micrOMEGAs 
2.0 code [20] with the MSSM model files modified to incorporate the mixed sneutrino. 
Superpartner and Higgs particle spectra are calculated using SuSpect [21]. We assume that 
only a single mixed sneutrino has a significant relic abundance today. Even if multiple 
right-handed sneutrinos are appreciably mixed with the active ones, this will still be true 
provided the light sneutrinos are not highly degenerate. 

In our calculations we fix the values of the MSSM parameters at the weak scale. We 
take the input parameters in the sneutrino sector to be the mixing angle 9, the LSP mass 
mp 1 , and the soft mass-squared for the left-handed sleptons, m\. Once 9 and are fixed, 
both the relic abundance and direct-detection rate (discussed in the following section) are 
both quite sensitive to m\. This is because increasing m 2 L increases the A (or X)-parameter, 
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thereby enhancing the annihilation rate via s-channel Higgs exchange and the cross-section 
for Higgs-mediated z/x-nucleon scattering. The gaugino masses Mi and M 2 can also be 
important in determining the relic abundance, as t-channel neutralino exchange is another 
potentially significant annihilation channel for v\. 

In Fig. [2] we display regions in m^-sinfl space that yield a relic abundance consistent 
with cosmological observations, for various values of Mi, M 2 , and m\. The other MSSM 
parameters are fixed as fi = 300 GeV, tan/? = 10, rriA = 500 GeV, mi R = (300 GeV) 2 , 
rriq = vn 2 UR = rn? dR = M 3 = (1 TeV) 2 , and A t — — 1 TeV, giving a Higgs mass of 116 GeV. 
Also shown are the constraints from the measurement of the invisible width of the Z and 
from the recent results from the XenonlO direct-detection experiment [8]. 

Let us consider the plots of Fig. [2j In plots (a)-(c), M\ and M 2 are held fixed as three 
different values of are used. Importantly, there are dramatic differences in the direct- 
detection constraints depending on whether elastic scattering via Z-exchange is suppressed 
- in which case Higgs exchange dominates - or unsuppressed. We will discuss the circum- 
stances in which the Z-exchange contribution is suppressed in the next section. 

If the Higgs-exchange contribution dominates, we see in (a)-(c) that v\ masses above the 
threshold for W + W" production and near the Z or Higgs poles are consistent with what we 
have learned about the dark matter abundance and with the latest XenonlO results. As 
is increased, the interesting regions in parameters space shift to smaller values of sin#. If 
instead the scattering via Z-exchange is unsuppressed, only the Higgs pole region is viable. 

In (d) we focus on a particularly light spectrum (both for sneutrinos and gauginos). In 
this case one sees that without the Z-exchange contribution, sneutrinos with the appropriate 
relic abundance are allowed over the entire mass range. With unsuppressed Z-exchange scat- 
tering, one is forced into the light mass range [vriy ^ 10 GeV). The precise mass below which 
this scenario is viable is not entirely certain, as the issue is sensitive to the highest velocity 
particles in the halo, for which a modified Gaussian is probably not a good description. 

In (e), we take tjil to be very large. As a consequence, there is a large A-term for the 
same value of sin#, making the s-channel Higgs annihilation more efficient, and allowing 
reasonable relic abundances for low values of sin (9. Consequently, a broad range of masses 
is viable, regardless of whether the scattering off of nuclei is dominated by Z- or Higgs- 
exchange. 

Finally, in (f), we consider the effect of modifying the width of the Higgs. In various 
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FIG. 2: Constraints on the sneutrino parameter space from requiring the correct relic abundance, 
from direct-detection experiments, and from the invisible Z-width measurement. The regions below 
the Z-width and direct-detection contours are allowed. The values taken for m\, Mi, and M2 are 
indicated in the plots, and the other MSSM parameters are as given in the text. 
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recent proposals [22j [23j |2U [251 12S1 |27J |28], the Higgs width is dominated by final states 
other than bb. This possibility is motivated by the fine tuning problem associated with 
raising the Higgs mass above the LEP limit. For our purposes, the importance of non- 
standard Higgs decays is that they can modify the form of the Higgs pole. We illustrate 
this by increasing the Yukawa coupling of the b-quark by a factor of five. One can see that 
this modification significantly impacts the allowed ranges for mixed-sneutrino dark matter 

- comparing (e) and (f), the mass ranges from 40 — 50 GeV and 60 — 80 GeV open up. 
It is worth emphasizing that the uncertainties regarding the decays of the Higgs can have 
significant consequences for the allowed parameter space of any dark matter model which 
involves annihilation through an s-channel Higgs. 

It is interesting to note that for the parameters used for Fig. |2](cl) , with ~ 10 GeV and 
sin 8 chosen to give the preferred dark matter abundance, the lightest Higgs boson decays 
invisibly, to v\i>*, more than 80% of the time. If is lowered further, the Higgs also 
develops an appreciable branching ratio into v-\p2 final states. Depending on the spectrum 

- for example, on whether decays to Xi are accessible - v 2 may decay dominantly to V\Z* . 
In this case, Higgs decays to v\V2 would most often produce a rather nondescript final 
state, although 20% of the time the Z* would decay to neutrinos, giving an additional 
contribution to the invisible width of the Higgs. In Fig. [3j we take Br{p2 v±Z*) = 100%, 
mi = 105 GeV, and mo 1 = 10 GeV , and plot the branching ratios for a 116 GeV Higgs to 
decay (i) invisibly and (ii) directly to standard model states, as functions of sin^. These 
branching ratios do not sum to unity because they exclude Higgs decays to z>i z> 2 with z> 2 
decaying visibly. In fact, we see that for the masses chosen for Fig. |3j and for the values of 
sin 9 preferred by cosmology, this third class of decays dominates. 

C. Direct and Indirect Detection of Mixed Sneutrinos 

Cryogenic detectors such as CDMS and liquid noble gas detectors such as XENON have 
made considerable strides since mixed sneutrinos were originally considered. In this section 
we consider the cross section for v\ - nucleon scattering, paying particular attention to the 
assumptions built in. As indicated in the previous section, we find, consistent with [H] , that 
broad parameter regions remain viable for mixed-sneutrino dark matter, and that broader 
regions open up in the lepton-number violating case. 
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FIG. 3: Branching ratios for Higgs decays directly to standard model states, and to invisible final 
states via either v\v\ or z>iz>2, for the parameters indicated. 

Including only the Z-exchange contribution, the cross section for v\ to scatter off of nuclei 

is 

a = \(A - Z) - (1 - 4 sin 2 9 W )Z] 2 sin 4 6, (4) 

where fi is the v\ - nucleus reduced mass. As discussed above, this cross section exceeds 
experimental limits for most parameters ranges, with exceptions at small ra^ , near the Higgs 
pole, or with heavy left-handed sneutrinos. 

Lepton number violation in the sneutrino mass matrix lifts mass degeneracy between 
the scalar and pseudo-scalar components of the lightest sneutrino, and scattering via Z- 
exchange occurs inelastically (i.e., through a transition from the scalar to the pseudoscalar 
or vice- versa) [HI ESI EQ]. As a consequence, particles with velocities below 

(3 mm = J— — (^^ + 5) (5) 
V 2M N E R \ fi J 1 ; 

are incapable of scattering. Here, M N and E R are the target nucleus mass and recoil energy, 

and 5 is the mass splitting between the scalar and pseudoscalar. Given that there are not 

expected to be any particles in the halo with velocities above the galactic escape velocity 1 , 



1 Or, in our reference frame, v esc + v ro t, where v ro t is the net total velocity from motion of the Earth about 
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by dialing S large (of order 100 keV), one can evade all direct detection constraints. 

In this case, scattering from Higgs exchange still constrains the theory [3]. In the decou- 
pling limit {m A 3> m h ), the cross section for Higgs-mediated i>i-nucleon scattering is 



where is the nucleon mass, guNN is the Higgs-nucleon coupling, and ghv x v x is the coupling 
of the light Higgs boson to the LSP sneutrino, 

ghvxvx = cos 2(3 sin 2 6 H — -= sin (3 sin 26. (7) 

v \J2 

The value of gtNN is subject to rather large uncertainties. Here we adopt the up and down 
quark, strange quark, and heavy quark contributions to this coupling given in Refs. [ST] , |32j, 
and [33J, respectively. This yields Qunn — 1-26 x 10~ 3 . Written in terms of this reference 
value, the cross section is 

i 



a = ( 9hNN - ) ( 9h ^ ) ( 115 GeV ) (2.48 x 10- 43 cm 2 ). 



(8) 



1. Relation to neutrino mass 

In the lepton-number- violating case gaugino loops generate neutrino mases [5j [6] . In the 
regime in which the diagram with a pure Wino running in the loop dominates, the correction 
to the neutrino mass is 

_ g 2 sin 2 ^m 1 m 2 v 
v ~ 32tt 2 ^ h31 [) 

where 

_ m 2 w mf \og[m 2 w /mf] + m 2 m 2 log[m 2 /m 2 ] + m\m) log[m 2 /m^] 
13 {™%-ml)(m\-m])(m 2 -m]) 
my are the masses of the light and heavy complex eigenstates, and 6 is the splitting between 
the scalar and pseudoscalar components of v\. 

For much of the parameter space, a splitting of order 100 keV is necessary to ensure 
that inelastic scattering at XENON is kinematically impossible. With the large mixings 
(and thus large A-terms) necessary to achieve the appropriate relic abundance, this mass 
splitting generates a neutrino mass of orderl eV. Combined with the limits on neutrino mass 



the sun, and of the sun about the galactic center. 
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from cosmology [31], which are roughly 1 eV for the sum of the neutrino mass this possibility 
seems excluded. This has been emphasized recently by [5]. In particular, the authors of [§] 
argue that inelasticity consistent with neutrino mass bounds change the allowed parameter 
space very little. Here we review some important caveats to this. 

One point is that the standard value used for galactic escape velocity, 650 km/s, may be 
too large. Indeed, the most recent simulations of Milky Way type galaxies [35] produce lower 
escape velocities (~ 450 km/s). Using the distributions of these simulations, rather than 
the much older halo parameters typically used to set limits, the allowed parameter regions 
for mixed-sneutrino dark matter do expand. For example, the parameters = 100 GeV, 
sin 9 = 0.18, mi — 300 GeV, M 2 = 400 GeV, and 5 = 50 keV give a realistic relic abundance 
and a direct-detection rate that is borderline at XENON, but only generate a neutrino mass 
of 0.32 eV. Although this mass pushes up against the cosmological limits, this example 
illustrates how the impact of 5 on direct detection rates is highly sensitive to assumptions 
about the halo. 

A second point is that the radiatively generated neutrino mass is suppressed as the 
neutralino masses are increased, whereas the v\ annihilation rate is insensitive to these 
masses if it is dominated by s-channel Higgs exchange. For example, for the same parameters 
as in the previous paragraph, raising Mi to 1 TeV does not change the relic abundance 
significantly, but does reduce the neutrino mass to 0.15 eV. 

Another possibility is that there might be an enhanced annihilation rate at smaller values 



of sin 9. For example, if the neutrino Yukawa couplings are large, as described in section [HA 



there are additional contributions to the annihilation rate coming from the v\V\h coupling 
2A 2 vcos 2 #. As this coupling can be parametrically comparable to A sin 29, the annihilation 
via s-channel Higgs can be considerably enhanced, even at smaller mixing angles. As pre- 
viously discussed, we leave the analysis of the model with large neutrino Yukawa couplings 
for future work. 

Finally, the neutrino mass of equation ^ requires a Majorana mass insertion. If the 
gauginos are Dirac, as described in [36], a radiative mass will not be generated. In these 
scenarios, one can disregard the radiative neutrino mass entirely, even for large S. 

When we consider the collider signatures of mixed sneutrinos, some of the parameter 
points we will study can accommodate mixed-sneutrino dark matter only if the scattering 
via Z exhange is strongly suppressed, possibly leading to a radiative neutrino mass that 
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is too large. Our main motivation for studying these parameter points is that LHC signa- 
tures for mixed sneutrinos are of interest in their own right, independent of the connection 
to dark matter. However, we also believe that because of the myriad astrophysical and 
particle physics uncertainties, a liberal take on which regions of parameter space may be 
cosmologically interesting is warranted. 



Indirect-detection experiments can also place important constraints on mixed-sneutrino 
dark matter; see [9] for a thorough discussion. Here we focus on the scenario in which 
inelasticity is relevant. 

If v\ particles are captured by the sun at a large enough rate, high-neutrinos produced 
in their decays can be detected on Earth. We will allow for the possibility that capture 
via Z exchange is suppressed by the mass splitting between the scalar and pseudoscalar 
components of v\, so that the capture rate is determined by the contribution from Higgs 
exchange. In this case we find, following [37], that the most stringent bounds from indirect 
detection experiments are not competitive with those from direct detection. For example, 
taking the parameters used for Figs. |2^b), with = 100 GeV and sin# chosen to give the 
desired relic abundance, we find a flux of upward through-going muons that is almost two 
orders of magnitude below the limits given in [38J. With the parameters used for Fig. [21(d) , 




and taking = 10 GeV, the sneutrinos now annihilate directly to neutrinos, but the 
predicted flux is still around an order of magnitude or more below current limits, depending 
on the flavors of neutrinos produced in the annihilations. 

As noted earlier, if the Higgs decays principally in non-standard fashion, then annihilation 
into neutrinos would be similarly suppressed if annihilation occurs through s-channel Higgs, 
although this certainly depends sensitively on the decay products of the Higgs boson. 



In supersymmetric theories in general, an intriguing possibility is that the true LSP is 
the gravitino, but that the lifetime of the NLSP is sufficiently long that the dark matter 
relic abundance is determined entirely by the freezeout of the NLSP [39, 40J. Typically, 



2. Indirect Constraints 





D. LSP gravitinos 
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the NLSP is imagined to be a stau or neutralino, but a mixed sneutrino could similarly 
serve as NLSP, decaying harmlessly into neutrino-gravitino. The sneutrino NLSP case was 
considered within the MSSM in [TP] . 

In this gravitino-LSP scenario the parameter space opens up dramatically, including 
regions with otherwise too large relic abundance, or regions where the XENON limits would 
have excluded mixed- sneutrino dark matter. This gives us extra motivation to be open- 
minded when studying the collider phenomenology of mixed sneutrinos. Unfortunately, this 
gravitino LSP scenario leads to no signals at dark matter detectors, either direct or indirect. 



III. LHC SIGNATURES 



What experimental signatures for mixed sneutrinos might be observed these at the LHC? 
This depends on the superpartner spectrum and, crucially, on the amount of mixing between 
the sterile and active sneutrinos. In section HI A we consider the case where the mixing 



angle 9 is quite small. In this case the mostly-sterile sneutrinos will be produced only 
rarely in cascade decays, unless they are the lightest superpartners. If they are the lightest 
superpartners, they will be produced in the decays of the NLSP, and the collider signatures 



depend on the identity of that particle. In section HI A we will see that if the NLSP is a right- 
handed slepton, a distinctive opposite-sign dilepton signature potentially emerges. If the 
NLSP is instead Xi> the collider phenomenology will be the same as with a neutralino LSP, 
but even in this case, there is a simple point to be made: a given cosmologically disfavored 
point in MSSM parameters space may become cosmologically viable with a mixed-sneutrino 
added at the bottom of the spectrum. 



In sections |IIIB| and |III C[ we consider additional signatures that become possible if the 
mixing angle is larger, 9 ^ 0.1. We have seen that, for these larger mixing angles, there is 
tension between having the correct relic abundance, evading direct-detection experiments, 
and satisfying the neutrino-mass bound. While in certain variations of the model this tension 
may be eliminated {e.g. in a scenario with Dirac gauginos), we regard these mixed-sneutrino 
signatures as important to study independent of whether the sneutrinos produced are the 
cold dark matter, for the following reasons: 

• Discovering mixed-sneutrinos at the LHC would shed important light on the nature 
of neutrino masses; it would suggest that the neutrinos are of Dirac type, or else, that 
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the seesaw scale is not be much larger than the weak scale. 

Even if the sneutrinos produced at the LHC are not the dark matter, they could still 
be relevant to the dark matter question. Here are two possible scenarios that illustrate 
this point: (1) The lightest mixed-sneutrino is the NLSP. It freezes out in the early 
universe and then decays to gravitino dark matter. (2) In addition to the sneutrinos 
with large enough mixing angles to be produced at the LHC, there is a lighter one 
with a smaller active component, suitable to be the cold dark matter. 



In section III B we consider kinematic edges that can appear in jet-lepton invariant mass 
distributions when charginos decay directly to mixed sneutrinos. We also discuss the pos- 
sibility of using mass estimation techniques to distinguish this signature from MSSM ones. 



In section HI C we study signatures associated with decays of the heavier sneutrinos to the 



lighter ones, involving Higgs and Z bosons 



A. Dilepton mass distributions 

Assuming that a mixed-sneutrino is the LSP, its presence at the end of every cascade 
decay chain makes for SUSY signals rich in leptons. In this section, we consider cascades 
involving a right-handed slepton NSLP, and find the following: 

• Opposite-sign same-flavor (OSSF) dileptons can arise from a two-body decay followed 
by a three-body decay, whereas in the MSSM they typically come from two-body 
followed by two-body, or from a single three-body decay. If the final state leptons 
are mostly fj, or e, this allows us to distinguish the mixed-sneutrino scenario from the 
MSSM over significant regions of parameter space. 

• If fi is produced, rather than e.R or a prominent signal is still possible and may 
be distinguishable from the MSSM. 

If the LSP is a weakly mixed sneutrino and the NLSP is a right-handed slepton, we have 
the following possible decays for the lightest neutralino: 

X \^~ Vl V X?-T!T X?-W, (11) 

where I — e, fi. The direct decay to V\V is suppressed by the small mixing angle, leaving 
f\T and IrI as the competing decay channels. We will assume for now that x? —* IrI has a 
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substantial branching ratio, and consider the case where Xi ~" > completely dominates at 
the end of this section. 

Under this assumption, we expect a large number of right-handed sleptons to be produced 
at the LHC. How do they decay? The possibilities are 

Ir -> v{W l R -> t x t\ l R -> z>xz//. (12) 

Even if the two-body decay to i)iW is kinematically accessible, it is not only suppressed by 
the sneutrino mixing angle, but also vanishes in the absence of left-right slepton mixing. 
For 6 ~ 0.05, it is typically negligible compared to the three-body decays for / = e, and 
only potentially competitive for I = /i. Note also that the relevant coupling here depends 
on flavor issues - if the active component of v x is entirely third-generation, these decays are 
absent. The second decay, to fir/, may or may not be kinematically allowed. Even if it is 
allowed, its kinematical suppression can easily make the third decay, to the dominant 
one. As we will illustrate by example below, this is true even though the decay to V\vl is 
mixing-suppressed. 

Assuming, then, that Xi ~^ W an d Ir — > V\vl both have substantial branching ratios, 
an opposite-sign same-flavor (OSSF) dilepton signature results. The / + /~ invariant mass 
distribution is predicted to have a kinematic endpoint at 

rnT+T- = m x\\[ l ~ KW^ 1 " ("W m fJ 2 - ( 13 ) 
In supersymmetric models, OSSF dilepton signatures, with associated kinematic edges, are 
quite common. How distinctive is the dilepton mass distribution in the case with a mixed- 
sneutrino LSP? 

In standard SUSY models, an OSSF dilepton signature can arise via the sequence of two- 
body decays x° — *■ — * (Xi^)^ + - One important difference compared to this standard 
case is simply that in the mixed-sneutrino LSP case, the dileptons come from Xi decays. 
Provided the relevant branching ratios are sizable we would thus typically expect a larger 
lepton multiplicity than in the case where X2 initiates the decays. However, Xi decays can 
also produce an OSSF dilepton signature in a scenario with a gravitino LSP and a right- 
handed slepton NLSP, through the sequence Xi ~^ (Ir)1 + ~^ (Gl~)l + (prompt decays of the 
NLSP slepton are possible for a low SUSY-breaking scale). 

These scenarios are easily distinguished from the mixed-sneutrino case by their dilepton 
invariant-mass distributions. The two-body/two-body sequences X2 ~^ (0^ + — * (Xil~)l + 
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and Xi ~~ * (Ir)1 + - * {Gl )l + both have the distribution 

dP 

oc mi+i-. 



(14) 



dmi+i- 

If we take the matrix element of the three-body decay to be constant and just consider the 
phase-space dependence, the two-body/three-body sequence xi ~^ (Ir)1 + ~~ y {pxvl^)^ gives 



dP 

dx 



oc x ( 1 — x 2 — fi 2 



In 



x 



(15) 



where we have defined x = mi+i-J \J m \o — m f , v — m l R l m xl^ anc ^ ^ = m v/ m i R - Normal- 
ized plots of these very different looking distributions are shown in Figure |4j 




FIG. 4: OSSF dilepton invariant-mass distributions from the sequence of two-body decays \% ~^ 

— ► (Xil~)l + (dashed), and from the two-body/three-body sequence x? ~~ ► ~~ ► 
(solid). For the latter, the amplitude of the three-body decay is set to be constant, and we take 



m? = 0.8 m^o and ms 



0.5 m^o. 

xi 



A softer dilepton invariant-mass distribution arises in the MSSM if the leptons come from 
a three-body decay such as x° — > Xil + l~ ■ in this case the endpoint of the distribution is 
just the mass difference between two neutralinos, 



m ™* =m x o-m x o. 



(16) 



Again taking the matrix element of the three-body decay to be constant and considering 
the phase-space dependence alone, the dilepton invariant-mass distribution is 

dP 



dx 



OC X 



^(l-x 2 )((l-K 2 ) 2 -x 2 ), 



(17) 



where x = mi+i- /(m^o — m^o) and K = 2m^o/(m^o — m^a). In the massless-LSP limit, this 



distribution and the two-body /three-body distribution of equation (15) both reduce to 

dP 



dx 



cx x(l 



x 



(18) 
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FIG. 5: OSSF dilepton invariant-mass distributions for the three-body decay xi ~* Xil + l~ 
(dashed), and for the two-body /three-body sequence Xi ~~ *■ (Ir)1 + ~^ (pivl~)l + (solid). For the 
latter, the amplitude of the three-body decay is set to be constant, and we take = 0.8 m^o. 
The distributions are identical for the massless-LSP case, shown in the middle. 

where x = mi+i- /m[+?f. Although these distributions are identical in the massless-LSP limit, 
they shift in opposite directions as the LSP mass increases, as shown in figure [5] There we 
see that the mixed-sneutrino distributions are significantly softer than what the MSSM 
three-body decays give, unless the LSP neutralino mass for the MSSM case is unusually 
small. Other observables, such as the overall lepton multiplicity, would likely help to further 
distinguish particular points in the parameter spaces of these two scenarios. Moreover, as 
we will see, the observed value of the kinematic endpoint of the mi+ l - distribution may not 
be easily reconciled with a very small value of m^o / m^o . 

We have used Monte Carlo simulation to explore the distinguishability of these scenarios 
further. We modified the SUSY-HIT package [H] for calculating superpartner masses and 
decay branching ratios to incorporate mixed sneutrinos and the associated decays. For 
example, the three-body decays of right-handed sleptons were implemented by appropriately 
modifying the the matrix elements given for three-body squark decays. We generated events 
with full SUSY production using Pythia 6.4 [12], with the pytbdy.f file modified to include 
the momentum dependence in the amplitudes for the three-body decays of right-handed 
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sleptons via off-shell bino, 
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\M(l R - 
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Xi 



m lu 
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mz — mf- 



(19) 
(20) 



Here mi v and m u c, are the lepton-neutrino and neutrino-sneutrino invariant masses, respec- 
tively. The momentum dependence for three-body neutralino decay amplitudes is already 
included in Pythia. After generating the fully showered and hadronized Pythia events, we 
passed these to the PGS 4.0 detector simulator [33], taking the granularity of the calorimeter 
grid to be 8r] x 8<p = 0.1 x 0.1. 

For the mixed-sneutrino case, we start with the mSUGRA-like high-scale parameters 
(10 GeV) 2 , Mi /2 = 450 GeV, tan/? = 10, A t = -500 GeV, and A b = A T = 0. Then 



T-,2 



m 



we add a weakly mixed sneutrino with 9 = 0.05 and mo 1 = 58 GeV to the bottom of the 
resulting spectrum. The physical superpartner masses for this point are given in Table [T| 
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TABLE I: Superpartner and Higgs boson masses for the parameter point used to study the dilepton 
signature in the mixed-sneutrino case. All masses are in GeV. 

The most important masses for the dilepton signature are m^o, m iK , and m^, whose 
values lead to a kinematic endpoint of m^L = 63 GeV. For the parameter point chosen, the 
branching ratios for the lightest neutralino are Br(xi — > IrI) = 41%, Br{\\ ~^ ^i r ) — 57%, 
and Br(xi — > V\v) = 2%. The branching ratios for the right-handed selectrons are are 
Br (en v\ve) = 97% and Br (en — > f\TV) = 3%, and the branching ratios for the right- 
handed smuons are Br(jl R — > UiUji) = 69%, Br(fl R — > UiW) = 29%, and Br(ft R — > f\TV) = 
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FIG. 6: Left: OSSF dilepton invariant-mass distributions for SUSY events and tt background. 
Right: flavor-subtracted opposite-sign dilepton invariant-mass distributions. 

2%. The total widths for and Jar are both hundreds of eV, so we can assume that they 
decay promply. Note that although we mix in three generations of sterile sneutrinos, all 
with the same mixing angle, the detectability of the dilepton signature does not rely on this 
simplifying assumption. For example, if there is only a single sterile sneutrino which mixes 
with the stau sneutrino alone, the branching ratios for jiR — ► v\V[i and e.R — > v\ve are both 
above 90% - the branching ratio actually goes up for p, R because it can no longer go to V\W . 

For this parameter point, we find using Prospino 2.0 [44J that the NLO cross section for 
squark and gluino production [451 HE] is 2.7 pb. We generate 80,000 events, corresponding to 
~ 30 fb _1 . In our analysis, we demand either an e + e" pair or a pair, where the leptons 

are required to have px > 10 GeV and \r)\ < 2.4. To suppress standard-model background, 
we further require YIiVt > 1500 GeV, where the sum is over jets with px > 20 GeV, leptons 
and photons with px > 10 GeV, and missing px- After this cut the leading standard- model 
background is from it, which we also simulate using Pythia. We generate 18.9 x 10 6 it 
events, corresponding to 23 fb" 1 of integrated luminosity, where we take a = 830 pb for the 
tt production cross section at NLO |47]. After the cuts, we have 850 it events after rescaling 
to 30 fb _1 , compared with 8,812 events from SUSY production. 

In the first plot of figure [6j we show the OSSF dilepton invariant-mass distribution for 
events passing the cuts. The distribution rapidly decreases as the mass approaches the 
expected endpoint mV+°l- = 63 GeV from below. It then levels off due to the relatively large 
SUSY background. This background can be dealt with using a standard flavor subtraction, 
as shown in the second plot of the same figure. To make that plot, the analysis is redone, this 
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FIG. 7: Flavor-subtracted OSSF dilepton invariant-mass distributions in (left) the gravitino LSP 
case and (right) the MSSM case with the three-body decay \% ~^ Xi^~- 

time demanding either an e + /i~ pair or a fi + e~ pair. The opposite-sign dilepton invariant 
mass distributions for these events are then subtracted from the OSSF distribution. We 
see that this procedure does a good job reducing the SUSY background, and the kinematic 
endpoint is evident, within a few GeV or so of the expected position at 63 GeV. 

A simple modification to this parameter point is to remove the mixed-sneutrino LSP with 
a light gravitino, with the assumption that e# and fiR decay promptly to it (this will not be 
the case if decays to fi are accessible). If we adjust the right-handed slepton mass to keep 
the OSSF dilepton endpoint near 63 GeV, we find the distribution shown in the first plot 
of figure [7j As expected, the differences compared to the the mixed-sneutrino case are clear 
due to the two-body kinematics of the relevant decays. 

We now want to compare the mixed-sneutrino distribution with that from an MSSM 
parameter point in which the three-body decays X2 Xil + l~ are important. So, we choose 
low-scale parameters such that the most relevant physical masses are m^± = 101 m^o = 101 
GeV, m^o = 37 GeV, mj L = 129 GeV, and = 102 GeV. The parameters are chosen so 
that the splitting between m^o and m^o gives a kinematic endpoint very close to the 63 GeV 
value from the mixed-sneutrino case. Also, the gaugino masses are taken as small as they 
can be consistent with negative results from direct SUSY searches at LEP. The rationale for 
doing this is to soften the dilepton invariant-mass distribution as much as possible, to see how 
similar to the mixed-sneturino distribution it can look (recall that the distributions should 
be very similar in the massless-LSP limit). Another way to soften the mi+i- distribution is 
to make the intermediate slepton in the three-body decay just barely off-shell. So, we make 
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FIG. 8: Flavor-subtracted OSSF dilepton invariant-mass distributions for SUSY events and tt 
background, in the mixed-sneutrino case. A second kinematic endpoint is visible at ~ 140 GeV. 

rriL small; however, we want to ensure that the two-body decay x® ~^ vv is not kinematically 
accessible, which prevents us from lowering arbitrarily close to m^o. Note that for the 
parameters chosen is just above m^o. 

The flavor-subtracted mi+i- distribution for this set of parameters is shown in the second 
plot of figure [7| The distribution is softer than in the gravitino LSP case, but still easily dis- 
tinguished from the mixed-sneutrino distribution of figure [6} Had the value of the kinematic 
endpoint for the mixed-sneutrino been significantly larger it would presumably allow one to 
find a point in the MSSM parameter space that gives a more similar looking distribution. 
For this reason we stress that other observables, such as the overall lepton multiplicity, may 
also be useful for distinguishing particular points in the two parameter spaces. A detailed 
study of a number of relevant observables at once would likely be efficient at ruling out 
candidate parameter points. 

For example, for the particular point we have chosen on the mixed-sneutrino side, an 
additional distinctive feature evident in the rrti+i- distribution is a second kinematic edge 
at ~ 140 GeV. This is shown in figure [HJ which is essentially a zoomed-in view of figure |6j 
going out to larger invariant masses. The flavor-subtracted distribution doesn't average to 
zero beyond ~ 63 GeV because OSSF dileptons can also come from the two-body/two-body 
sequence x° ~" > — * (x?^ ~K + > an d the kinematic endpoint for this sequence is indeed 
near 140 GeV for the parameters chosen. The possible presence of two edges, one associated 
with X2 decays and the other associated %i decays, is one more puzzle piece one could use 
to distinguish points in mixed-sneutrino parameter space from points in MSSM parameter 
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space. 

Decays to r's 

We conclude this section by considering the possibility that Xi ~ > ^iT is the dominant 
decay of the lightest neutralino, with nearly 100% branching ratio. In this case, the signa- 
ture involving e + e~ or pairs is absent, and the situation becomes more challenging 
experimentally. If fi decays dominantly to V\VT, then one could hope to observe an endpoint 
in the ditau invariant mass distribution. However, if fi — > v^W is kinematically accessible, 
it is likely to dominate. Here we assume this two-body decay is kinematically accessible and 
consider the detectability of Xi — ¥ iji) T+ — > (^i^)t + , with the W decaying leptonically. 

To study this issue we take the high-scale parameters fh 2 = (50 GeV) 2 , My2 = 350 GeV, 
tan/? = 10, A t = —500 GeV, and Aj, = A T = 0. Then we add a weakly mixed sneutrino 
with 9 = 0.05 and mo 1 = 51 GeV. For the purposes of this study, the most important 
masses are m^o = 141 GeV, = 134 GeV, and mo 1 = 51 GeV, and the most important 
branching ratios are Br{x\ — > f\r) = 94% and Br(?\ — > v{W") = 94%. In the sequence 
Xi —> (^i)t + - * (^i[H^ _ ])r + — > (b>i[lv])T + , the invariant-mass of the final-state lepton and 
tau has an upper bound of 37 GeV for the masses considered. 

We find the total SUSY production cross section for this point to be 10.6 pb for pp 
collisions at y/s = 14 TeV, and work with ~ 587,000 events, corresponding to about 55 
fb _1 of integrated luminosity. We take all events with at least one reconstructed r and at 
least one isolated lepton (e or /x), and as before we put apply a cut J2pt > 1500 GeV to 
reduce standard-model backgrounds. The invariant mass m T i is calculated for all r-lepton 
pairings with opposite sign, giving the OS m T i distribution shown in the first plot of figure 
[9| Repeating the same procedure, this time requiring the r and the lepton to have the 
same sign, gives the SS distribution in the same plot. The peak at low invariant-mass in 
the OS distribution arises from events with Xi ~ > (^i) r+ ~~ > {vi[W~])t + — > (vi[Iv])t + . The 
OS distribution is similar to the SS distribution beyond this peak, but an excess in the OS 
distribution over the SS distribution does persist beyond the expected 37 GeV endpoint, 
because there are other ways to produce opposite-sign leptons and taus in association with 
each other —* (^i) r+ ~" > (^i[W~])r + — > (ui[W})t + being just one example). 

To gauge whether this signature would be observable above background, we include the 
tt sample generated with Pythia. Taking the difference between the OS and SS distributions 
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FIG. 9: Left: opposite-sign (OS) and same-sign (SS) lepton-tau invariant-mass distributions from 
SUSY production alone. Right: opposite-sign minus same-sign distributions for SUSY plus ti 
production, and for ti production alone. 

for both the SUSY and ti backgrounds, we find the subtracted distributions shown in the 
second plot of figure [9} Despite the significant ti and SUSY backgrounds, a rather dramatic 
fall-off in the distribution is evident for invariant masses between 30 and 40 GeV. 

In the MSSM, the decay v T — > ?iW + may occur due to the left-right stau mixing. If the 
stau subsequently decays as T\ — > x°t, then we have associated W-r production, just as 
we had in the mixed-sneutrino case. Here we do not explore in detail the extent to which 
this MSSM decay sequence could be distinguished from the mixed-sneutrino decay sequence 
considered above, but simply note that in the mixed-sneutrino scenario the signal has the 
potential to be much more prominent, given that it originates from Xi decays rather than 
v T decays. 

B. Jet-lepton mass distributions 

At the LHC, most SUSY evens will begin with squark or gluino production, and so 
leptons produced in association with mixed sneutrinos will be typically be accompanied by 
hard jets. In this section, we consider the impact of mixed sneutrinos on jet-lepton invariant 
mass distributions, and find the following: 

• The decay chain q — > (x ± ) ( l ~^ leads to a prominant edge in this distribution, 

providing a potentially distinctive signature for mixed-sneutrino production. 



350 - 

300 - 

250 - 

£L 200 - 



50 - r 



25 



• A recently proposed mass-estimation method [H] can be used to probe the spectrum 
of the theory when applied to events involving this decay chain. That method can 
thus be used to help distinguish mixed-sneutrino and MSSM scenarios, although in 
our implementation it does not reliably estimate the sneutrino mass. 

If the sneutrino mixing angle 9 is large enough, then chargino and neutralino decays 
directly to v\ can become important. In fact, it is easy to imagine a situation in which 
x\ v\v and xt ~^ both have nearly 100% branching ratios. To illustrate this with 
a concrete example, we take the high-scale parameters fh 2 = (200 GeV) 2 , m# Hd = 
Mi/a = 300 GeV, tan (3 = 10, A t = -500 GeV, and A b = A T = 0. We add to the resulting 
MSSM spectrum mixed sneutrinos with 9 = 0.2 and = 108 GeV. The superpartner 
spectrum for these parameters is given in table III} From this table we see that the only two- 
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TABLE II: Superpartner and Higgs boson masses for the parameter point used to study the jet- 
lepton signature. All masses are in GeV. 

body decays available to xf are to ul, XiW, and T\v. The branching ratios for these final 
states are 95%, 4%, and < 1%, respectively. As before, we make the simplifying assumption 
that all three generations of sterile sneutrinos have equal mixing angles with the active 
states. If there is only one sterile sneutrino, which only mixes appreciably with z> r , then 
X will decay almost exclusively to u%t. In this case detecting the chargino decays becomes 
more challenging. 

Although v\ is lighter than Xi m the above spectrum, this ordering is not important for 
the signature we are about to explore. If we had taken = 140 GeV instead, we would 
still have Br(xt &l) — 93%. So it is not essential for the following discussion that V\ 
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is the LSP. We should also note that, given how dominant Xi ~^ &l is for the parameters 
we've chosen, that decay mode still easily dominates for more modest values of the mixing 
angle, 9 ~ 0.1. So, signatures associated with xf ~^ can exist for those smaller mixing 
angles as well. 

For the parameters chosen, x® a ^ so decays dominantly straight to ui, with a branching 
ratio greater than 99%. Because of the presence of ui, both x\ an d X2 appear as missing 
energy in the detector. In particular, although there are leptons produced in chargino decays, 
there is no OSSF dilepton signal initiated by x\ decays. 

Chargino production via squark decay leads to the following sequence: q — > (x^Q — > 
(pil)q. One then expects to observe a kinematic endpoint in the jet-lepton invariant- mass 
distribution at 

m™r = rmy/i - (m xi /m,)^l - (m^/m^)*. (21) 

When we simulate events that include this sequence of decays using Pythia, the chargino is 
decayed isotropically in its rest frame, and angular correlations between the lepton and jet 
are lost. Taking xf t° be pure charged wino, which it nearly is for the chosen parameters, 
the squarks that can initiate this sequence are ul, cIl, u* l , or d* L . If these are produced 
with equal abundance, then the quark-lepton angular correlations average out to zero, even 
if we focus on a particular sign for the charge of the lepton. In this case the fact that 
Pythia does not keep track of the angular correlations is not important, and the true m q i 
distribution looks the same as the distribution for a 2-body/2-body sequence of decays with 
an intermediate scalar, shown as the dashed line of figure |4} 

However, there is no reason to expect that u L , d^, u* L , and d* L will be produced in 
equal abundance. For example, for the sample point chosen above, the dominant SUSY 
production is squark +gluino. The gluino decays with roughly equal probabilities to all four 
of these possibilities, but the parton distribution functions dictate that the squark produced 
is less likely to be di than ul, and less likely still to be an anti-squark. The sequence 
Ul —> {x + )d —>■ (vil + )d gives a quark and a lepton with opposite helicities, while the sequence 
di — > (x~)u —> {v*l~)u gives a quark and lepton with the same helicity. So, the jet and 
lepton tend to be more back-to-back when produced by d^, and more in the same direction 
when produced by ul- Given that ul is produced more abundantly than di, we should then 
expect that the combined m q i distribution will look somewhat softer than the dashed line of 
figure [|J without the sharp edge. On the other hand, if we focus on jet - l~ invariant-masses, 
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then we should expect that distribution to be even harder than for the case without angular 
correlations. This is because l~ is produced in the sequences beginning with (which gives 
a quark and a lepton with the same helicity) or with u* L (which gives an antiquark and a 
lepton with the opposite helicity). Assuming di is produced more abundantly than u* L , the 
jets and leptons will then tend to be more back-to-back on average. 

Although these angular-correlation issues are important, we set them aside in what fol- 
lows. The jet-lepton kinematic edge we identify below may be softened when both signs 
of lepton charge are allowed, but by requiring negatively charged leptons, an even harder 
distribution should result. 

For the parameter point chosen above, we find a total SUSY production cross section of 
20.4 pb. We generate ~ 160,000 events, corresponding to ~ 8 fb _1 , and keep events with 
the following characteristics: 

• Exactly two jets with px > 150 GeV. 

• Exactly one isolated lepton with pt > 10 GeV. 

• A transverse mass my > 250 GeV. 



Missing transverse energy 1$t > 250 GeV. 
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FIG. 10: Left: invariant-mass distribution for jets and leptons in the same event, and rescaled 
invariant-mass distribution for jets and leptons in different events. Right: the subtracted distribu- 
tion. 



We find that the number of tt events passing these cuts is more than a factor of 20 
smaller than that from SUSY production. Using Alpgen 2.12 [35], we estimate the VF+jets 



28 



background by obtaining W+2 jet events with a generator-level cut on the two jets, pr > 100 
GeV. We find that this source of background is even more suppressed than the tt background. 



In figure 10 we show the jet-lepton invariant-mass distribution for the SUSY events, 
where for each event we include the two invariant masses obtained by pairing the isolated 
lepton with both of the hard jets. A fairly steep drop-off in this distribution is seen near 
expected kinematic endpoint, which is at 568 GeV. There is a large combinatorial background 
associated with incorrect lepton-jet pairings, whose shape one can attempt to guess by 
looking at the invariant-mass distribution for leptons-jet pairings from different events. A 
similar method is used to study SUSY ditau signatures in [SO], for example. Whether 
the distribution for incorrect lepton-jet pairings from the same event and the distribution 
for lepton-jet pairings from different events are similar should obviously depend on the 
strength of the correlations in the momenta of incorrectly paired jets and leptons from 
the same event. Because the squarks will tend to be at least somewhat back-to-back, one 
would expect correlations at some level. However, in the present example, the shape of 
the distribution for lepton-jet pairings from different events matches rather well with the 
shape of the same-event distribution beyond the expected kinematic endpoint. Rescaling the 
different-events distribution to match the same-event distribution at high invariant mass, 
and then subtracting that rescaled distribution off, the drop-off near the kinematic endpoint 



becomes clearer, as shown in the second plot of figure 10 . The bump in the distribution at 
relatively low invariant mass (rriji < 300 GeV) is due to cascades involving t\. 

A jet-lepton invariant-mass distribution of the sort shown can also arise without mixing, 
from production of ordinary sneutrinos. Furthermore, a jet-lepton signal can be produced 
in other ways in the MSSM. A first example, relevant if the decay \\ ~~ * x\W~ is dominant, 
is given by the sequence q — ► (Xi)q — * (XiW~)q ) with W~ decaying leptonically. A second 
example is given by the decay sequence q — * (xi)q ~- > (^l]* 7 )*? — * (\xV]^)Q- I n both of these 
examples, the lepton and jet are accompanied by a neutrino, and so the kinematics are 
different than in case of sneutrino production. One would expect the jet-lepton invariant- 
mass distributions to reflect these differences at some level. In fact, for the second example 
one can show that the distribution is softer than for a 2-body/2-body sequence with an 
intermediate scalar, regardless of the sign of the lepton considered. The analysis for the 
first example is more complex, as there is no intermediate scalar in the decay chain, but 
here there are other things to go on as well. For example, one could look for events where 
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the decay Xi ~^ XiW occurs on both sides of the event, with one W decaying leptonically 
and the other hadronically. If a significant number of hadronic W's were reconstructed by 
looking for events of this type, the xf ~^ Xi^ ± interpretation would have to be favored 
over the x ~^ &il interpretation. Another potentially important difference is that one has 
flavor universality for the decays W~ —>■ W, but not necessarily for the decays x ~~ * 

Other observables may be useful for distinguishing production of mixed-sneutrinos and 
ordinary sneutrinos. For example, for the parameter point considered above, x% decays 
invisibly, and the flavor-subtracted OSSF dilepton invariant-mass distribution has no par- 
ticularly distinctive features. On the MSSM side, for much of the parameter space for which 
Xi — > v*l occurs, x\ a ^ so occurs, giving rise to an OSSF dilepton signature upon the 

subsequent decay l* L — > These decays tend to come along with each other because in 

the MSSM the masses of v and II are split only by electroweak symmetry breaking, 

m\ -ml<m 2 w . (22) 

Provided that m^± is not much larger than m^o, the decay X2 ~^ ^1} thus tends to be 
kinematically accessible when Xi ~ * v*l is. The essential point is that, because the mass of 
the mixed sneutrino is not directly linked to the mass of charged slepton, it is easier than 
in the MSSM to have signals for sneutrino production in the absence of signals for charged 
slepton production. 

One way to have Xi ~^ v*l without X2 ~> ^1} m t fle MSSM is to have a closely-spaced 
spectrum with m^, < m-± and mj L > m^o. Mass measurements would clearly be helpful 
in distinguishing a mixed-sneutrino scenario scenario from this MSSM one. For example, 
if it were established that the sneutrino-chargino mass splitting were quite large, it would 
disfavor the MSSM scenario just described. 

Mass Estimation 

Here we consider how the technique proposed by [H] can be used to probe the mass 
spectrum. The authors of that paper consider a general situation in which the sequence 
Y — > IX and X — > I'N occurs on both sides of an event. The particle N is invisible, so the 
final state topology involves four leptons and missing energy. A typical SUSY example of 
this situation has Y = X21 X = I, and iV = Xi- F° r the mixed-sneutrino scenario considered 
in this section, we have a large number of events with q — > x ±c l ~^ &lq on either side. The 
event topology is thus quite similar, but with two of the four leptons replaced with quarks. 
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For a set of candidate values for the unknown masses my, mx, m^, one can check whether 
the observed kinematics of a given event of this type are consistent with those values. The 
procedure of [IB] is to fix two of the masses and keep track of the number of allowed events 
as the third mass is scanned. A candidate value for the third mass is identified by looking 
for a a dramatic feature in the in the resulting distribution, e.g. a sharp drop or peak in the 
number of allowed events. In our implementation we apply a smoothing procedure to the 
events distribution and identify the candidate mass as the point where the second derivative 
of the distribution is minimized. 

These steps are then iterated - the third mass is fixed at its new candidate value while 
the first mass is scanned, and so on. It was found in [48J that this procedure does not 
converge, and it was suggested that the actual masses can be estimated as the ones that 
give a global peak in the number of consistent events as the iterations are performed. In our 
implementation we find that for some of the events samples the candidate masses quickly 
settle near a final value, while for others the candidate masses continue to jump around 
indefinitely. Even in the case where the candidate masses continue to jump around, they at 
least wind up in stable "orbits" after a sufficient number of iterations. We take the average 
values of these orbits as the mass estimates for a given event sample. 

Although we do not find that our implementation of this procedure leads to a reliable 
estimate of the sneutrino mass, it does give a reliable estimate of the chargino-sneutrino mass 
splitting. We first apply the technique to the same parameter point described above, taking 
the most optimistic case where the chargino decays produce electrons and muons and not 
taus. We select events with exactly two jets with px > 150 GeV, and exactly two leptons 
with pt > 10 GeV. The leptons are required to have the same sign in order to suppress 
standard model backgrounds. We find that the efficiencies for passing these cuts are 1.4% 
for the SUSY sample and 1.3 x 10~ 5 for the tt sample, giving a ratio of SUSY events to it 
events of over 20. More problematic than the tt background is the background from SUSY 
events in which the selected jets and leptons do not come from the desired decay chains. 
Using 12 sets of 1,000 events, we obtain 12 estimates for (m^, m~±, m^). Combined, these 
event samples correspond to an integrated luminosity of roughly 42 fb _1 . For each event 
sample, we take the average values of the candidate masses at large iteration number as the 
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mass estimates for that sample. Averaging these estimates gives 

(m 9 ~,m^±,m^) = (688 ± 33 GeV, 239 ± 27 GeV, 110 ± 30 GeV) 
m~± - m Pl = 129 ± 7 GeV, 

compared with the actual values, (m^, m^±, m^) = (684 — 688 GeV, 229 GeV, 108 GeV) 
and m~± — rap l = 121 GeV. If the chargino decays produce e, //, and r with the same 
probability, the resolution worsens somewhat. Performing the same procedure on 10 sets of 
1,000 events, we obtain 

(m q -,m-±,m^) = (720 ± 40 GeV, 263 ± 35 GeV, 140 ± 42 GeV) 
m .± - mH = 123 ± 10 GeV. 

If we now redo the same analysis with the same parameters except with the sneutrino mass 
increased to 142 GeV, we obtain mp 1 = 116 ± 22 GeV and m^± — = 92 ± 7 GeV 
(versus an actual splitting of 87 GeV). Increasing the sneutrino mass further to 185 GeV 
gives = 124 ± 23 GeV and m-± — = 48 ± 2 GeV (versus an actual splitting of 44 
GeV). 

These results show that although the sneutrino mass estimates do not follow the actual 
values closely, the estimates of the chargino-sneutrino mass splitting do. So, this analysis 
can be used to find evidence against a closely-spaced spectrum of the type described above. 
A more sophisticated implementation of this mass-estimation method, or different methods 
such as those proposed in [SH E2] may do a better job at estimating the sneutrino mass 
itself. 

C. Associated production of Z/h with a lepton 

For some regions of parameter space, the heavier sneutrino would be produced at the 
LHC. In particular, if a chargino with significant charged-wino component is heavier than 
the heavier sneutrino, it can decay to that sneutrino and a charged lepton. The question 
is then how the heavier sneutrino decays. One possibility is v<i — > x\ u i which is just what 
one might expect in the MSSM. Here we consider a more distinctive scenario, in which the 
splitting with the lighter sneutrino and the sneutrino mixing angle are both large enough 
that the branching ratios for z>2 — > viZ and z>2 — > V\h are both significant. In this case we 
find 
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• Cascade production of Z bosons leads to a clean signature if the sneutrino contains 
e or \x flavor. The shape of the trilepton invariant mass distribution can in principle 
distinguish this signature from MSSM signatures involving Z bosons. 

• If the sneutrino is dominantly r flavored, a much larger integrated luminosity of the 
LHC is necessary to see cascade Z production 

• Cascade production of Higgs bosons can also lead to distinctive signatures, both 
through the h — > bb and h — > 77 decay channels. 

As an example of a point in parameter space with these interesting heavy-sneutrino 



decays, consider the weak-scale values shown in table III, and the resulting mass spectrum 



of table IV We find that the SUSY production cross-section for this spectrum is 24 pb, and 
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TABLE IV: Superpartner and Higgs boson masses for the parameter point of table III 



the relevant branching ratios for the signals that interest us here are Br(xt ~ > ^0 = 32% 
(I = e, /i), Br{y2 — > V\Z) = 37%, and Br{yi — > v\K) = 37%. These decays lead to a 
trilepton signature for the case of Z production, and Ibb and £77 signatures for the case of 
Higgs production. 
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To explore these signatures, we generate 640k SUSY events, corresponding to about 27 
fb _1 of integrated luminosity. We select events with the following properties to study the 
lepton-Z signature: 

• Three isolated leptons with > 10 GeV. Two of these must be of opposite sign and 
same flavor, with \mi+i- — mz\ < 10 GeV. 

• Y2pt > 800 GeV, where the sum is over jets with px > 20 GeV, leptons and photons 
with p T > 10 GeV, and missing p T . 

These cuts select 1,323 SUSY events. We find the background from ti and WZ production 
(also generated using Pythia) are negligible by comparison. The trilepton invariant mass 



distribution for events passing the cuts is shown in the first plot of figure 11 The distri 
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FIG. 11: Left: invariant-mass distribution for trileptons in the same event, and rescaled invariant- 
mass distribution for Z candidates and leptons in different events. Right: the subtracted distribu- 
tion. 



bution shows a steep drop-off near the expected endpoint at 364 GeV. Along the lines of 
what was done for the jet-lepton mass distribution, one can attempt to subtract off the 
SUSY background by guessing that the shape of the SUSY background contribution to the 
distribution should be similar to to the distribution obtained by matching the Z candidates 
from one event with leptons from different events. Subtracting off this distribution gives the 



distribution shown in the second plot of figure [TTJ A lower kinematic endpoint should in 
principle be observed at around 148 GeV. With greater statistics it is possible that a second 
edge near this value would become clear in the subtracted distribution. 
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Z bosons can also be produced in association with leptons in the MSSM, in processes 
such as X2 — > {Xi)Z — > (p*l)Z. In principle, this MSSM process is distinguishable from 
the mixed-sneutrino signature considered above by examining the Z-l invariant-mass dis- 
tribution. In the mixed-sneutrino case, the intermediate 1)2 decays isotropically in its rest 
frame. The amplitude-squared for the chargino decay into i/\lZ is thus constant, and the 
mzi distribution is trapezoidal, rising linearly between lower and upper kinematic edges. In 
the MSSM sequence xi ~^ {Xi)Z — > (v*l)Z, there is no intermediate scalar, and there will 
generically be some dependence of the amplitude-squared on mzi- The chargino-chargino-Z 
coupling can be written as 

X2Y(clPl + c r P r )xi, (23) 

where the couplings cl and cr are determined by the mixing in the chargino sector. These 
couplings are typically not equal. In the narrow-width approximation, the amplitude- 
squared for the decay X2 v*lZ takes the form a + /3(|cl| 2 — \cR\ 2 )mf z , where a and 
/3 are constants. Provided the chargino- Z coupling is indeed chiral, the invariant mass dis- 
tributions thus differ in the MSSM and mixed-sneutrino cases. A similar issue was raised in 
considering the possible production of Z bosons in top-squark decays in [53]. 

If only the v T mixes with a sterile neutrino, the signature becomes much less clean than 
what we have considered. In this case v e ^ decay straight to %\v, and so one is forced to look 
for signals from associated r — Z production. To explore these we use a much larger sample 
of SUSY events corresponding to roughly 170 fb _1 , and impose the same cuts as before, 
except that his time we require a reconstructed r and two opposite sign, same flavor leptons 
that reconstruct a Z . After these cuts, we are left with 3,020 SUSY events. We estimate 
that the it background gives fewer events by a factor of 4.8. In our plots we will not include 
the it contribution, as we have not generated a large enough it sample; it is likely that more 
carefully chosen cuts can improve the quoted signal to background ratio. 

In the first plot of figure [12] we show the tII invariant-mass distribution for events passing 
the cuts, along with the (rescaled) distribution obtained by matching Z candidates from 
one event with r's from a different event. The second plot of the same figure shows the 
subtracted distribution, which does have something that looks like an endpoint. The excess 
of events in the different-event distribution at low invariant mass is presumably due to the 
fact that the isolation criteria for r's and leptons in the same event for that are not enforced 
for r's and leptons in different events. 
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FIG. 12: Left: rll invariant-mass distribution for t's and Z candidates in the same event, and 
rescaled invariant-mass distribution for r's and Z candidates in different events. Right: the sub- 
tracted distribution. 

To study the lepton-Higgs signature, we first select events with the following character- 
istics: 

• Two 6-tagged jets with p T > 20 GeV. 

• Exactly one isolated lepton with p T > 10 GeV. 

• A transverse mass my > 200 GeV. 

• Missing transverse energy J^t > 200 GeV. 

• TIvt > 800 GeV, where the sum is over jets with px > 20 GeV, leptons and photons 
with pt > 10 GeV, and missing p?. 

Working with the same 27 fb _1 sample used to analyze the lepton-Z signature, a total of 2,432 
SUSY events survive these cuts, a factor of eight larger than the number of ti events that 
survive. The invariant mass distribution for events passing these cuts is shown in the first 



plot of figure 13 A peak, although not a clean one, is evident at around 100 GeV, somewhat 
below the actual Higgs mass of 114 GeV. We keep events events with |m&& — 100 GeV| < 20 
GeV, and the m&y invariant-mass distribution for those events is shown in the second plot 



of figure 13 The distribution shows a significant drop-off near the expected endpoint, at 
369 GeV. 

For the parameter point chosen, the rate of Higgs production is large enough that it the 
Higgs can also be seen through its decays to photons. For this analysis, we impose the same 
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FIG. 13: Left: invariant-mass distribution for pairs of 6-tagged jets. Right: for events in the 
peak of the m^b distribution, invariant-mass distribution obtained by pairing Higgs candidates and 
leptons. 

cuts as for the bbl analysis, except that we require two photons with pt > 10 GeV instead 
of two 6-tagged jets. A total of 60 SUSY events pass these cuts. We have checked that 
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FIG. 14: Left: invariant-mass distribution for pairs of photons. Right: for events in the peak of the 
m 77 distribution, invariant-mass distribution obtained by pairing Higgs candidates and leptons. 



the backgrounds from ti production and W / 77+jets (estimated using Alpgen) are negligible. 
(In fact, it is possible that it would be advantageous to relax the kinematic cuts.) The 77 
invariant mass distribution for events passing these cuts is shown in the first plot of figure 
14 , and has an extremely clear peak. The m 77 z distribution for events in this peak is shown 



in the second plot of figure 14; it falls below ~ 400 GeV as one would expect, but more 



statistics would be required to learn much from this distribution. 
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IV. CONCLUSIONS 

In this paper we have explored the cosmology and collider phenomenology of mixed 
sneutrinos. Recent progress made by direct-detection experiments makes mixed-sneutrino 
dark matter quite constrained, but still viable. In the absence of lepton-number violation, 
parameter regions that give an appropriate relic abundance while evading direct-detection 
constraints include (1) light sneutrinos, mo 1 < 10 GeV, with relatively light gauginos and 
a rather large mixing angle, 9 ~ 0.3, (2) small mixing angles, 9 < 0.07, with near the 
Higgs funnel, or (3) small mixing angles with above threshold for annihilation to W 
pairs, for large values of heavier sneutrino masses, m f2 > 500 GeV. 

Lepton-number violation in the sneutrino mass-sqaured matrix can suppress scattering 
of v\ off of nuclei via Z exchange, making somewhat larger mixing angles viable. How- 
ever, this lepton-number violation produces radiative contributions to neutrino mass. If the 
lepton-number violation is large enough to dramatically suppress the elastic scattering via Z 
exchange, the radiatively generated mass tends to approach or exceed the upper bound from 
cosmology. In principle one can suppress this radiatively generated neutrino mass by making 
the gauginos heavy, while still achieving a realistic relic abundance through annihilations 
mediated by Higgs and Z exchange, or by making the gauginos Dirac. 

We have studied LHC signatures of mixed sneutrinos in general, without requiring the 
mixed sneutrino responsible for the signal to be the dark matter. If the mixed-sneutrino is 
the LSP, with a very small mixing angle, then the NLSP will will typically be the only particle 
that decays to it with a large branching ratio. If this particle is the lightest neutralino, then 
the only effect of the mixed sneutrino is to alter the connection between collider physics 
and cosmology. If the NLSP is instead a right-handed slepton one expects an unusually 
large lepton multiplicity in the SUSY signal. Moreover, decays of the lightest neutralino 
can lead to an interesting opposite-sign, same-flavor (OSSF) dilepton signature. Since the 
dilepton signature arises from a two-body decay followed by a three-body decay, the shape 
of the dilepton invariant-mass distribution is significantly different from the sequence of 
two-body decays \\ — * — * (Xi^)^ + - It can also be distinguished from that arising 
from a three-body decay ~^ Xi^~ \ depending on factors such as the observed kinematic 
endpoint of the distribution. If the lightest neutralino decays to f\ rather than cr or jin, the 
experimental signatures become more difficult to extract, but there is still the possibility 
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of observing a large excess of events with taus produced in association with opposite-sign 
leptons, where the leptons come from W bosons produced in the decay f x — > v{W . 

For larger mixing angles, it is possible that the mixed sneutrinos will be produced copi- 
ously at the LHC through chargino decays, or even from the decays of heavier sneutrinos. 
In the first case, the sequence q — > (x^Q ~ * {p\l)q gives rise to a kinematic edge in the jet- 
lepton invariant-mass distribution. For broad regions of parameter space, not only — > x ± 
but also X2 decays dominantly to z>i, and in these regions the jet-lepton signature is present 
in the absence of an OSSF dilepton signature. This situation can also arise due to ordinary 
sneutrino production in the MSSM, but we have shown that mass-estimation methods may 
be helpful for distinguishing the mixed-sneutrino and MSSM scenarios, due to the fact that 
the sneutrino - charged slepton mass splitting is an electroweak symmetry breaking effect in 
the MSSM, but not in the mixed-sneutrino case. 

Finally, if the predominantly left-handed sneutrinos z> 2 are produced at the LHC, the 
decays z> 2 v\Z and u 2 — > v\h may be important if kinematically accessible. Events 
with z> 2 is produced from chargino decay may then have Higgs or Z bosons produced in 
association with leptons. In this case, a distinctive bbl invariant mass distribution can arise 
from Higgs production, and a distinctive trilepton invariant-mass distribution can arise from 
Z production. Because these decay chains feature an intermediate scalar, z> 2 , they can in 
principle be distinguished from MSSM decay chains such as X2 ~^ (Xi)Z ~^ (v*l)Z, where 
the chargino-chargino-Z coupling is in general chiral. We leave a detailed analysis of this 
issue for future work. 
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